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According to the Gottesmann-Knill theorem the non-Gaussian states are necessary component for 
a nontrivial quantum computation. We show two efficient and deterministic methods of x'''' non- 
Gaussian state generation for a cavity mode using a single trapped ion. Both require ion motional 
state transfer to the cavity field. The methods are experimentally feasible. The first is based on 
the well-known protocol for an ion finite motional superposition state generation. It allows for an 
arbitrary good approximation of x''^^ non-Gaussian states. We give criteria based on the Wigner 
function which quantify the error resulting from the approximation. The second and novel method 
enables an exact non-Gaussian state generation using one laser pulse only. 



I. INTRODUCTION 

The Gottesmann-Knill theorem for continuous vari- 
ables (CV) states that quantum computing based only 
on components described by one- or two-mode quadratic 
Hamiltonians, Gaussian states input, and measurements 
on canonical variables can be efficiently simulated by a 
classical computation. That means, although some of the 
algorithms are of fundamentally quantum nature, they 
do not provide any speedup over classical processes 
It has also been proved that construction of a CV uni- 
versal quantum computer for transformations that are 
polynomial in those variables requires cubic or higher 
non-linear operations Ql- Therefore, investigation of the 
non-Gaussian transformations and states generation is 
crucial for a nontrivial quantum computation. More re- 
cently it has been noted that Kerr-like nonlinearities, in a 
variety of systems, enable high precision quantum mctrol- 
ogythat would otherwise require entanglement to achieve 
it % 

Over the last decade ion trap experiments have led 
the emerging technologies of coherent quantum control, 
especially in quantum information theory [J] and quan- 
tum computation Q . Those applications require efficient 
creation and precise manipulation of both electronic and 
motional trapped ion state. Various theoretical propos- 
als on how to produce nonclassical arbitrary states of ion 
motion have been discussed. In experiment, Fock number 
states [1] , coherent states [? [, vacuum squeezed states [1] , 
and Schrodinger cat states 3| have been realized. In this 
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latter case, the state is an entangled state of the vibra- 
tional and electronic degrees of freedom. On the contrary, 
in this paper we give a deterministic way to prepare the 
vibrational degree of freedom in a non-Gaussian state 
that is not entangled with the electronic states. Accord- 
ing to our knowledge, neither non-Gaussian state (other 
than an entangled cat state) nor a superposition of more 
than two coherent states has been observed so far. 

In the case of photons there is no practical method of 
non-Gaussian state generation so far. Efforts have been 
made to explore a class of the x^^'' non-Gaussian states 
produced using a photon coherent state \a) interacting 
with x'-^'' Kerr nonlincarity in an optical fiber 



|vl/(a,r))=e-^^-^ 



e'2 



(1) 



This class of non-Gaussian states, parametrized by the 
unitless evolution parameter r [lO|, is the most popu- 
lar one. The state ([1]) is known also as the Kerr state. 
In general, this is a highly nonclassical state and after 
a certain time of evolution r in the fiber its Wigner 
function would take negative values in the phase space 
However, the nonlincarity in a fiber, or any other 
experimentally achievable Kerr medium, is too small. 



^{3) ^ 10- 
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, to reach a highly nonlinear regime and 
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thus produce the negativity in an experimentally reason- 
able time, before it is destroyed by dissipation [l^]. Al- 
though microstructured fibers seem to be more promising 
with x^'^'' — 10^^^ their length does not exceed Im 
with current technology. 

The most known examples of the Kerr state are the cat 
states e~'^^/^\ia) -\- e*'^/'^| — ia) corresponding to r = tt, 
which have been found useful for studies of quantum de- 
coherence and quantum-classical boundary [3, 13 • The 



2 



larger cat states for which the two components \ia) and 
I — ia) are ahiiost orthogonal {a > 1.5) find their appli- 
cation in quantum information processing [l5l | and quan- 
tum computation p^ . [l7j . 

Recently, there has been introduced a probabilistic 
method of non-Gaussian state generation relying on a 
conditiona.1 photon subtraction from a squeezed vacuum 
state [H, [13, HO, mi- Such a state is a good approxi- 
mation for the cat state if the amplitude is small a < 1. 
However, it is neither practical to subtract more than two 
photons in the experiment nor to produce a state with 
a > 1.6 [2l[. Moreover, the state is not pure. 

Another approach is suggested by recent achievements 
in the the Polzik group ^22] based on their demonstration 
of teleportation of a quantum state from optical to matter 
degrees of freedom. If this process was reversed, so that 
a non-Gaussian state is either teleported or mapped from 
matter to light, optical Kerr states might be robustly gen- 
erated and thus overcome obstacles which are met with 
in trying to generate optical non-Gaussian states directly. 
Of course this requires us to demonstrate a way to make 
non-Gaussian states in matter degrees of freedom. The 
proposals for transferring an arbitrary motion al q uantum 
state of an atom to a cavity field already exist |23l. [2^ . [25l . 
The effort to integrate ion traps and optical fields [2y| 
might offer a path based on the extraordinary level of 
coherent quantum control one has over the vibrational 
degrees of freedom for trapped ions. 

In this paper we discuss two efficient and determin- 
istic methods of x^^^ non-Gaussian state generation for 
light using a single trapped ion. Both methods require 
ion motional state transfer to the cavity mode. The first 
method is based on the well-known protocol [27j for an 
ion finite motional superposition state generation. It al- 
lows to produce the x^'^^ non-Gaussian states with ar- 
bitrary good approximation. We give criteria, based on 
Wigner function comparison and its measurement preci- 
sion, which quantify the error resulting from the approx- 
imation. The second method is novel and it enables an 
exact non-Gaussian state generation using one laser pulse 
only. We point out that a Wigner function measurement 
of ion motional state can be performed using currently 
available technology and already existing experimental 
schemas. We also suggest a quantum metrology applica- 
tion, based on the work of Caves and co workers [3]. 

This paper is organized as follows. In section |TT] we 
show that, applying a well-known protocol, one can pro- 
duce an approximated x*-'^^ non-Gaussian state of motion 
for a trapped ion. The method relies on using a series 
of laser pulses to couple electronic and vibrational de- 
grees of freedom to effect the desired state preparation 
for the motional degree of freedom. We give the criteria 
for quantifying the extent to which the prepared state ap- 
proximates the desired non-Gaussian state, and discuss 
the technical limitations of the method. In section Hill we 
present an alternative, and novel, method for generation 
of an exact x*-^' non-Gaussian state. We also discuss the 
range of application of this method. We finish the pa- 



per with conclusions and a brief discussion of possible 
applications to quantum metrology. 



II. SERIES OF LASER PULSES METHOD 



An ion in a Paul trap [28 
bitrary state of the form 



may be prepared in an ar- 



\^)^6\^'^)\e)+f3\^a)\g), (2) 
using the method proposed in [23, , where 
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are finite superpositions of ion motional states {M < 00), 
\n) is a Fock state of a harmonic oscillator potential 
in the trap, \g) and |e) are the ion electronic ground 
and excited states respectively. Parameters S and /3 are 
complex numbers obeying the normalization constraint 
|5p -I- = 1 and are set at the start of the experiment. 

The method is based on applying a series of alternating 
laser pulses tuned, first to the carrier, and then the red 
sideband, transition of the trapped ion. It works both 
within and beyond the Lamb-Dicke regime. The ion is 
initially prepared in its electronic and vibronic (motional) 
ground state \0,g). Adjusting the Rabi frequencies and 
duration of each pulse properly, one could achieve = 



Wn equal to 



1 



(4) 



fc=0 fe! 



These coefficients correspond to the coefficients of a 
quantum non-Gaussian state resulting from the unitary 
evolution of a self-Kerr interaction ^ decomposed in the 
Fock basis, up to the normalization factor (which results 
from the fact that we cut off the infinite sum in ^ and 
take into account only the first M terms). 

This method produces an approximated x^^'' non- 
Gaussian state |vI'(*^)(a,T)) corresponding to the state 
reached via a Kerr interaction with an arbitrary value of 
evolution parameter r. Furthermore the effective value 
of T can be much larger than can be achieved via unitary 
interaction under a realistic optical Kerr interaction were 
the bosonic degree of freedom an optical field mode. A 
proper choice of the cut-off value M enables one to ap- 
proximate the state ^ arbitrarily closely. If we set (5 = 
for simplicity, the ion state is given by 



|vl'(*^)(a,r)) 



1 n 

y ^e'5"("-i)|n)|5). 



(5) 



The number of required pulses for preparing the state ([5]) 
is equal to 2Af 

l^^^'^a^T)) =RmCm-i ■ ... •i?iCo|0,g), (6) 
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FIG. 1: The coefficients io„ of a x'^' non-Gaussian state de- 
composition in a Fock basis evaluated for a = 2 - the top 
figure and a — 5 - the bottom figure. The dots above the 
red dashed line take values greater or equal to 10"'^. li a — 2 
{a = 5) the significant coefficients range from wo (wb) to wia 

{W5l). 



where Cj and Rj denote a carrier and a red sideband 
laser pulse respectively. Therefore M should be as small 
as possible [30| . 

The hint for existence of a good non-Gaussian state 
approximation comes from a simple observation that only 
a finite amount of coefhcients Wn contribute to the sum 
P]) significantly and the exact number of them is strongly 
a dependent. The coefhcients, evaluated for exemplary 
values of the amplitude, a = 2 and a — 5, arc depicted 
on Fig. [TJ 

For a given value of an amplitude a the choice of the 
cut-off value AI is based on three criteria. All criteria rely 
on comparing the Wigner function W^'^-'^t = 27r,7,7*) 
of an approximated state ([5]) with the untruncated one 
W{t = 27r,7,7*) evaluated for ([T]). Please note that the 
comparison is made for r = 27r. We chose this particular 
value of the evolution parameter because for this value 
the original Wigner function is given by a simple analytic 
formula, which can be computed directly. 

The first criterion derives from an investigation of the 
approximated Wigner function isolines. The smaller the 
value of A/ the more the function deviates from an ideal 



Gaussian function for r = 27r: the isolines are no longer 
circles. The isolines of interest can be chosen arbitrarily. 
We however chose O.I, 0.3, 0.5 and calculate, for each 
isoline, the ratio between the most and the least distant 
point with respect to the point (a, 0) separately. For an 
ideal circle the ratio is always equal to I. 

In order to take into account the interference re- 
sulting from the approximate state, which takes the 
values around zero, we calculate the maximal and 
average ratio of the difference | |1^(t = 27r, 7, 7*) — 
W^'^'^^T = 27r,7,7*)| to the value of the ideal Wigner 
function, for all points 7 in a phase space. Therefore, the 
second criterion gives the maximal and average relative 
error of the approximation respectively. 

The third criterion reveals the percentage of points 7 
in a phase space for which II^(^^) (r = 2n, 7,7*) = W{t = 
27r, 7, 7*) for a given precision. We assume a precision of 
order of either ±10~^ or ±10~^ to be good enough, since 
it relates to the accuracy of both the Wigner function 
computer visualisation using the density plots and the 
Wigner function reconstruction using quantum tomogra- 

phy mi. 

The above criteria turn out to be more subtle than the 
well-known fidelity J- which measures the overlap of the 
approximated and the original state in the phase space 



M 

E- 

m=0 



_(7) 

Obviously, the fidelity approaches unity in the limit of 
M 00. 

We would like to point out that having the ion already 
prepared in state ([5]) one could measure the Wigner func- 
tion of its vibronic state u sing either the standard method 
of quantum tomography 13X1331 or the method of direct 
measurement developed in [34|. This would allow for de- 
tailed investigation of x'"^^ non-Gaussian states, which 
has never been verified experimentally so far. 

The direct Wigner function measurement method is es- 
pecially interesting. It relies on the fact that the Wigner 
function of a displaced vibronic state of an ion is related 
to the probability of finding the ion in the ground and the 
excited state. The probabilities are measured by detect- 
ing a fluorescence signal. Such a measurement has been 
done for a light in a cavity with accuracy of ±0.2 [ssj . 



Example: I*'"' (a = 2, r)) generation 

Estimation of the minimal number of laser pulses for 
|\I>(*-'^) (a = 2,r)) generation, which approximates the 
original state \'9{a = 2, r)) for any value of the evolution 
parameter r well, requires comparison of the approxi- 
mated Wigner functions W^'^\t = 27r,7,7*) evaluated 
for a few different cut-off values M, with the original 
Wigner function W{t = 271,7,7*). The selection of the 
best possible cut-off values is based on analysis of sig- 
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nificant coefficients Wn- For a = 2 only the coefficients 
ranging from wq to Wiq are greater than 10"'^. The co- 
efficients wio to u>i3 arc of order of 10^^. Therefore, we 
test 10 < M < 16. 

We next investigate the isohnes. Table U shows the 
ratios of the most to the least distant point, with respect 
to the point of (2,0) for the isolincs (0.1, 0.3, 0.5) of 
the approximated Wigner fmictions W^^^"* (r = 27r, 7, 7*). 
The minimal value of M for which the ratios evaluated 
for W'-'^^'>{t = 2tt, 7, 7*) are equal to the ratios evaluated 
for the original Gauss function using the grid (discretized 
phase space) with a step equal to A7 = 0.04 is equal to 
14. It means that for the given grid and greater values 
of M, the numerical simulations of the Wigner function 
will not differ. Therefore, the value of M = 14 could be 
regarded as an appropriate cut-off value. 

The series of 28 laser pulses is experimentally feasible 
jsot . However, decreasing the quality of approximation 
only a little, one can diminish the number of pulses to 20. 
The ratios evaluated for M = 10 (A/ = 9) and isolines 
0.1, 0.3 and 0.5 differ from unity of 20% (33%), 10% 
(23%) and 8% (19%) respectively. This is the minimal 
cut-off value for which the average error equal to is below 
1% (see table For M = 9 the error is equal to 

1.67%. The maximal error is equal to 6.39% (11.30%). 
Within a given precision 10~^ for M = 10 {M = 9) there 
are around 76 % (64 %) points for which the numerically 
obtained Wigner function values are equal to the values of 
the ideal function. For the precision of 10"'^ and M = 10 
(M = 9) there are 36 % (26 %) of such points. 

These results show that M — 10 is still a good ap- 
proximation for the Wigner function analysis and mea- 
surement. The approximated functions for Af = 10 and 
Af = 14 with marked isolines are depicted on Fig. [51 We 
also include the plot for M = 9 ior comparison. 

The fidelities evaluated for |«'(9)(a = 2,t)), |*(io)(a = 
2,r)) and |*(")(q; = 2,r)) are equal to = 0.9838, 

= 0.9943, and = 0.9999 respectively. 

Below we present the time of pulses duration and their 
phase required to generate an exemplary compass state 
|\I>(io)(q; = 2,t = 7i'/2)), which is a superposition of four 
coherent states 
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These quantities were evaluated for the following Rabi 
frequencies ilc = IMHz, fii? = 100 kHz corresponding 
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FIG. 2: The Wigner function for an approximated y^^'' non- 
Gaussian state: |^'^^-'(a = 2, r = 27r)) - the top figure, 
\¥^'^\a = 2,T = 27r)) - tlie top middle figure, \¥^^\ol = 
2, T = 2-k)) ~ the bottom middle figure, |*(a = 2, r = 2-k)) - 
the bottom figure. 
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FIG. 3: The Wigner function for an approximated x'"^^ non- 
Gaussian state: [^''■^^(a = 2,t = ^)) - the top figure, 
1*'^") (a = 2, T = f )> - the middle figure, |*(q = 2, r = 2ir)) 
- the bottom figure. 



to carrier transition and red sideband respectively, and 
Lamb-Dicke parameter r] = 0.02. The total time of state 
build-up is equal to tc — 6.7 ms. 

The Wigner function of the approximated compass 
state for M = 10 is is presented on Fig. [31 We also 
present the Wigner functions for M = 9 and the original 
one for |^'(a = 2, r = ^)) for comparison. 



Technical limitations 

At the end of this section we estimate the maximum 
value of the amplitude a for which the method works. 
The maximum value seems to be a ~ 2.3. 

This limitation does not result from the dccohcrcnce 
of the ion state. Both the time of single pulse duration 
and the total duration of state built-up remain within the 
coherence of the ion for the vibronic state 190 ms. The 
coherence time for the electronic state is equal to 1.4 ms. 

The number of laser pulses is limited due to the finite 
binding energy in the trap. The cut-off value M corre- 
sponds to the maximal ion excitation \M). The \M = 17) 
is the upper limit for trapping the ion so far [36| . 



III. ONE LASER PULSE METHOD 

The one laser pulse method of x^^-* non-Gaussian state 
generation requires one carrier resonance pulse applied 
to the ion cooled down to a Lamb-Dicke regime. The vi- 
bronic state of the ion is initially prepared in a coherent 
state \a). The advantage of this method over the first one 
is that it allows for generation of the original Kerr state 
([T|) without any approximations. However, not all val- 
ues of the evolution parameter are accessible: the pulse 
duration is limited due to the ion decoherence. 

The interaction between an ion and a laser pulse of 
Rabi frequency is governed by the following Hamilto- 



(8) 



where is an electronic state rising operator, a is a 
vibronic state annihilation operator, is a trapping fre- 
quency, and 77 is a Lamb-Dicke parameter. Using the 
expansion of the exponens function into a Taylor series 



^M:(ae— + ate-)^ (9) 
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in the rotating wave approximation and in an interaction 
picture the Hamiltonian is approximated by the first five 
terms of the expansion 
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The terms in the first line in the above formula govern 
the free ion evolution. The last term in (fT0|l we asso- 
ciate with a self-Kerr interaction Hamiltonian, up to the 
electronic state operators. Therefore, the unitary evo- 
lution operator resulting from the nonlinear part of the 
Hamiltonian is given by 



U{t) 
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It allows for reading out the effective value of the evolu- 
tion parameter 



(12) 



In the further discussion we neglect the electronic state 
evolution since it may be made to take no part in the 
dynamics. To do this one first needs to prepare the elec- 
tronic state in an eigenstate of using a 7r/2 pulse tuned 
to the carrier transition, but the pulse must not excite 
the vibrational degree of freedom in any way. This can 
be done by making the pulse propagate orthogonal to the 
vibrational axis. 

Please note that by way of contrast to the method 
described in the previous section this method is inde- 
pendent of the initial amplitude a. The laser pulse 
duration required for the Kerr state generation is the 
same for all values of the amplitude. In other words, 
the time required for generation e.g. |^'(q; = 2, r)) and 
\^{a = 5,t)) is the same. However, accessibility of the 
evolution parameter Tcff is limited instead. The greater 
Rabi frequency O and Lamb-Dicke parameter value are 
the shorter pulse duration t must be to obtain Toh. On 
the other hand, the ion has to remain within the Lamb- 
Dicke regime in order to hold the expansion (|10[1 true. 

Figure d] shows the time of the pulse duration assum- 
ing its Rabi frequency = 10 MHz and the Lamb-Dicke 
parameters 7] = 0.1, r/ = 0.3, r/ = 0.02, and ?] = 0.03. 
For the first two Lamb-Dicke parameters it is less than 
800 /iS, which is within the coherence time of vibronic 
ion state (190 ms), to obtain Toff ~ tt and therefore pro- 
duce the cat state. For the two remaining 77 values the 
duration of the pulse is of order of tenth of second. 

The duration of the pulse required for |5'(q!, r = 0.04)) 
generation is equal to t = 0.16 ms for — 0.1 and t ~ 
1.98 fis for T] = 0.3 (t = 0.1s for r] = 0.02 and t = 0.2 /xs 
for T] = 0.03). For r = 0.04 the first negativities in the 
Wigner function become visible [Tl] |. 

This method allows for obtaining also the other coher- 
ent state superpositions easily. For example, for Toff = ^ 
one achieves superposition of six states and if t^s = f one 
achieves superposition of four states (the compass state) 
[37I I. For 77 = 0.3 and |^'(q;, t = 7r/3)) the pulse duration 
is equal to t ~ 51.71 /is, and for |^'(Q;,r = 7r/2)) it is 
equal to t ~ 77.57 /is. These values are experimentally 
feasible. 

The one laser pulse method of x^^-* non-Gaussian state 
generation is limited due to introduced cut-off in the ex- 
pansion It is valid only if we take the expansion 
for a small parameter. To estimate the small parameter 
we approximate the quantum operators by the classical 
amplitudes a ~ a 



E 

A:=0 



kl 



2 a 



(13) 



and read out the small parameter in the expansion to be 
277a. The above series is decreasing if 770; < ^• 



This constraint is enough to ensure that the contribu- 
tion from the higher order terms is negligible. For exam- 
ple, for a = 5 and 77 = 0.1, obeying = ^ (?? = 0.09 and 

a < ^) the coefficients are equal to: 0.5 (0.405) 

for fc = 2, 0.041 (0.027) for fc = 4, 0.0013 (0.0007) for 
fc = 6. For coefficients obeying ct < ^ the sixth order 
coefficient is of two orders of magnitude smaller than the 
forth order coefficient. 



IV. WEAK FORCE DETECTION 

We now show how the cat-state weak force detection 
protocol proposed in Munro et al (ssj can be implemented 
in an ion trap given the ability to engineer an approxi- 
mate cat state as we have described. It must be admitted 
that there is no compelling case to use the vibrational mo- 
tion of a trapped ion as a weak force detector. However, 
the method we propose here would be a nice demonstra- 
tion of how non-Gaussian states can beat the standard 
quantum limit for weak force detection. 

Suppose one was able to prepare the vibrational degree 
of freedom in the non-Gaussian state 
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with the amplitude a real. When a weak force acts it can 
be described by the action of the unitary displacement 
operator D{ie) = exp(7ea^— iea) acting on the initial state 
to give the output state \il)o) = D{ie)\il]i). Using the 
result that 



giIm(ia£)|Q,. 
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D{ie)\a) 

we find that for e <C Q', and a real that 

\tlJo) = cos(7r/4 + ae)\+) +i sin(7r/4 + ae)\-) (16) 
where the even and odd parity states are given by 
1 



± 



V2 
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In other words, the weak force is well approximated by 
a rotation in the two dimensional parity subspace. As in 
[3^ it then follows that the minimum detectable force is 
then given by 
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2^ 



(18) 



which beats the standard quantum limit by a factor of 

Ion trap provide a simple way to reach this lower 
bound. Suppose that after the weak force has acted, the 
total vibrational and electronic state is given by Eq. p6p . 
The first step is to make a 7r/2 rotation of the electronic 
state, followed by the conditional rotation 



R 



(19) 
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which can easily be done [39|. Finally another 7r/2 rota- 
tion gives the state 

I^Ao) = cos(7r/4 + ae)|+)|g) +isin(7r/4 + ae)|-)|e). (20) 

The electronic state can now be readout, and the proba- 
bility to find the ion in, say, the excited state is 

P+ = i (1 - sin(2a6)) . (21) 

Sampling this distribution gives the estimation for the 
force with a minimum detectable force that is inversely 
proportional to a. For example, the minimum force re- 
quired to shift the interference distribution by one fringe 
is emm = 7r/(4a). 

V. CONCLUSIONS 

In the paper we have presented two experimentally fea- 
sible methods of x*-^-* non-Gaussian states generation for 
light using a trapped ion. The first method is based on a 
well-known protocol (a series of laser pulses) and allows 
for generation of an arbitrarily well approximated non- 
Gaussian state. The approximation is quantified by three 



criteria based on Wigner function analysis. This method 
is only limited by technical parameters of the trap: the 
binding energy. 

The second method enables an exact non-Gaussian 
state generation using one laser pulse. It is limited by the 
decoherence time of the motional ion state. However, ad- 
justing the laser pulse and trap parameters properly one 
is able to produce a cat state. 

Based on the proposed protocols [l^, [l^, [1^ we believe 
that non-Gaussian state transfer from ion motion to a 
light beam is possible in a foreseeable future. It will 
enable application of a nontrivial quantum computing 
protocols based for example on coherent states. This is 
also a step towards exploring so far unknown branch of 
quantum optics such as non-Gaussian states. 
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Isoline 


M 


9 


10 


11 


12 


13 


14 


15 


16 


0.1 




1.33 (33%) 


1.20 (20%) 


1.10 (10%) 


1.04 (-) 


1.01 (-) 


1.01 (-) 


1.00 (-) 


1.00 (-) 


0.3 




1.23 (23%) 


1.10 (10%) 


1.04 (-) 


1.01 (-) 


1.00 (-) 


1.00 (-) 


1.01 (-) 


1.01 (-) 


0.5 




1.19 (19%) 


1.08 (8%) 


1.01 (-) 


1.04 (-) 


1.04 (-) 


1.04 (-) 


1.04 (-) 


1.04 (-) 



TABLE I: The ratios of the most to the least distant point, with respect to tlie point (2,0), for given isolines 0.1, 0.3, 0.5, of 
the approximated Wigner function W^'^'\t = 27r,7,7*) for 9 < M < 16. 



Accuracy^^ 


M 


9 


10 


11 


12 


13 


14 


15 


16 










10-2 




64% 


76% 


90% 


99% 


100% 


100% 


100% 


100% 


10"^ 




26% 


36% 


44% 


53% 


65% 


80% 


97,% 


100% 



TABLE IL The number of points 7 in the phase space for which W^'^'\t = 27r,7,7*) — W{t = 271,7,7*) s-t a given accuracy 
for 9 < M < 16. 



M 

Erroi^^^ 


9 


10 


11 


12 


13 


14 


15 


16 


Average 
Maximal 


I. 67% 

II. 30% 


0.98 % 
6.39 % 


0.55 % 
3.48 % 


0.30 % 
1.84% 


0.15% 
0.94 % 


0.08 % 
0.47% 


0.04 % 
0.22 % 


0.02 % 
0.11% 



TABLE III: The average and the maximal error estimated for a — 2 and for 9 < A/ < 16. 
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0.15 



0.5 1 1.5 2 2.5 3 



FIG. 4: The time of pulse duration required for Kerr state 
|*(a,r)) generation for Lamb-Dicke parameter rj = 0.1 - the 
top figure, rj = 0.3 - tire top middle figure, 77 — 0.02 - the 
bottom middle figure, r; = 0.03 - the bottom figure. The red 
horizontal line denotes time duration of 1 s. 



